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A theory of resonant Raman scattering spectroscopy of one dimensional electronic systems is 
developed on the assumptions that (i) the excitations of the one dimensional electronic system are 
described by the Luttinger Liquid model, (ii) Raman processes involve virtual excitations from a 
filled valence band to an empty state of the one dimensional electronic system and (iii) excitonic 
interactions between the valence and conduction bands may be neglected. Closed form analytic 
expressions are obtained for the Raman scattering cross sections, and are evaluated analytically and 
numerically for scattering in the polarized channel, revealing a "double-peak" structure with the 
lower peak involving multispinon excitations with total spin S = and the higher peak being the 
conventional plasmon. A key feature of our results is a nontrivial power law dependence, involving 
the Luttinger Liquid exponents, of the dependence of the Raman cross sections on the difference of 
the laser frequency from resonance. We find that near resonance the calculated ratio of intensity 
in the lower energy feature to the intensity in the higher energy feature saturates at a value of the 
order of unity (times a factor of the ratio of the velocities of the two modes). We explicate the 
differences between the 'Luttinger liquid' and 'Fermi liquid' calculations of RRS spectra and argue 
that excitonic effects, neglected in all treatments so far, are essential for explaining the intensity 
ratios observed in quantum wires. We also discuss other Luttinger liquid features which may be 
observed in future RRS experiments. 

I. INTRODUCTION 

One-dimensional (ID) physics has become increasingly important to condensed matter science in recent decades 
because many systems, including for example single wall carbon nanotubesi, organic conductors^, superconducting 
narowires^, spin chains^, and even ultracold atoms in highly anisotropic magneto-optical traps^, have been developed 
to levels that allow experimental studies of hitherto unprecedented extent and precision. Among these ID (or quasi- 
1D) systems, the semiconductor quantum wire structure (QWR) is one of the most important and widely studied, 
because of its simple band structure and highly tunable doping density. Also, great progress in micro-fabrication 
techniques has made high quality samples available 7 . 

Theoretically, one dimensional systems are very well understood. The pioneering work of Tomonaga™, Luttinger— 
and Haldaneifi along with many other studies has produced an essentially complete understanding of the low energy 
physics. However, the relation between the theoretical results and experimental data is not as clear and direct as one 
would like. For example, only a few unambiguous observations of the fundamental concept of spin-charge separation 
have been reportedii*i2iiii. 

Resonant Raman Spectroscopy (RRS) has become a powerful tool for studying the elementary excitations of 
electrons in many different systems. Applications to low dimensional doped semiconductor nanostructures (two- 
dimensional quantum wells or one-dimensional quantum wires) have been particularly prominent 6,14 . In the usual 
RRS experiment, schematically represented in Fig. da), external photons are absorbed at one frequency and mo- 
mentum by exciting an electron from the valence band to the conduction band (step 1) and then emitted at another 
frequency and momentum via the recombination of the hole in the valence band with another electron in the conduc- 
tion band (step 2). The final state of the system contains one or more particle-hole pairs excited in the conduction 
band. The dispersion of these particle-hole pair states may be inferred from the energy and momentum difference 
between the incident and the scattered photons (Stokes shift), while the nature (collective mode (plasmon) and multi- 
pair excitation) may be inferred from the dependence of the scattering amplitude on the polarization of the incoming 
and outgoing photons and on the energy of the incident photon. 

The standard theory of Raman spectroscopy in semiconductorsi^i 6 ^*^, is based on the Fermi liquid (FL) quasi- 
particle picture and neglects resonance effects. The neglect of resonance effects means that in the unpolarized geometry 
(identical polarization of incoming and outgoing light) the Raman process couples to the electron density operator 
so the Raman cross section is proportional to the dynamical structure facto»i2iiS of the conduction band electrons. 
Within Fermi liquid theory and in the long wavelength limit relevant to Raman scattering, the electron structure 
factor has a delta function peak at the plasmon (CDE) energy and a much weaker feature at a lower energy associated 
with incoherent particle-hole pairs. The lower energy feature is referred to as the single particle excitation (SPE) 
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peak. The calculated RRS intensity therefore has strong spectral peaks when the energy difference between the 
incident and outgoing photons coincides with the collective CDE mode frequencies at the wavevector defined by the 
experimental geometry, while the intensity at the SPE energy is much weaker than that at the CDE energy (the 
ratio is of the order of the square of the ratio of momentum transferred by the light to the Fermi momentum of 
the electron gas ; for typical experimental relevant parameters, about three orders of magnitude). This theoretical 
result, however, is in qualitative disagreement with the experimental data2fi*2L22, in which the polarized spectrum 
exhibits a "double peak" structure with comparable-intensity peaks at both CDE and SPE modes. This puzzling 
feature22*2ii22 of an ubiquitous strong SPE peak in addition to the expected CDE peak occurs in one, two, and 
even in three dimensional doped semiconductor nanostructures, for both intrasubband and intersubband excitations. 
Many theoretical proposals22i24i2£i2§i2L2& have been made to explain this two-peak RRS puzzle. However, two of us 
have recently argued 2 ^ that within the standard theory, which neglects resonance effects, none of the above proposed 
modified mechanisms theory can even qualitatively explain the experimental data. 

But the situation changes when resonance effects are included. In an important paper, Sassetti and Kramer (SK)^4 
first proposed that in ID systems, the prominent lower energy SPE peak is due to " spinon" excitations whose coupling 
to light is enhanced by resonance effects. The qualitative idea that resonance effects can strongly affect the relative 
absorption cross-sections of different modes is indeed important. Unfortunately, as we have recently noted, the theory 
presented in Ref. 24] suffers from two technical flaws. First, it is not self-consistent: it uses a Fermi-liquid-based 
expression to account for the resonance effects, but a Luttinger-liquid-based expression to account for the conduction 
band dynamics. As we shall show below (and have already mentioned in a previous brief communication 2 ^) Luttinger 
liquid physics affects the matrix element in a crucial way. Second, while the SK calculation correctly notes that as 
resonance is approached the coupling becomes long ranged in space, it omits the equally important fact that the 
coupling becomes long ranged also in time. 

In this paper, which amplifies and extends our previous short communication 2 ^, we derive and present a complete, 
closed-form expression for the RRS scattering amplitudes in the Luttinger liquid model and calculate the resulting RRS 
spectra in different resonance conditions. We treat both the analytically tractable case of short ranged interactions 
and the physically relevant case of the Coulomb Luttinger liquid. In the Coulomb case we predict an asymmetric 
broadening of the spectral peak in the higher energy side arising from the curvature of the plasmon dispersion. We 
find that most aspects of the RRS spectra are similar to those predicted by the Fermi liquid approach. As noted in 
our previous work 2 ^, characteristic Luttinger effects are revealed in the dependence of the intensities on the difference 
of the laser frequency from the resonance condition. Going beyond the bosonic expansion developed in our earlier 
work 2 ^, we present explicit results for the total spectral weights of the charge boson and spin boson excitations 
in the polarized RRS channel. Far from resonance we find that the spin-singlet mode at energy, u> = qvf, has 
spectral weight much smaller than that of the charge boson (i.e. plasmon); however as resonance is approached the 
weights in the two contributions become comparable. We explain the difference between the Luttinger liquid results 
presented here and the 'Fermi Liquid' results presented previousljiiSiiSiiLi2iS2i^&. Our results remain inconsistent with 
present experimental data, so we argue that Luttinger liquid effects have not yet been unambiguously detected in 
RRS experiments. One possibility is that excitonic effects, neglected in the present and previous treatments, are 
important. 

The paper is organized as follows: In Sec. ^we develop a general theory of Resonant Raman Scattering in Luttinger 
liquids. We then apply our theory to calculate the spectrum for short- and long-ranged electron-electron interaction 
by using bosonic expansion method in Sec. IIIII In Sec. IIVI we go beyond the bosonic expansion and calculate the full 
spectral weights of the charge boson and spin-singlet excitations in the polarized channel. We compare our calculation 
to the previous 'Fermi liquid' calculations in Sec. E] We critically discuss the assumptions and resonance effects of 
our results and compare them with the present experimental data in Sec. I VII Finally we summarize our results in 
Sec. EH 

II. RESONANT RAMAN SCATTERING CROSS SECTION OF A LUTTINGER LIQUID 

In this section we present a derivation of the RRS scattering cross section for a Luttinger liquid. We consider an 
idealized model of a quasi-one-dimensional system, in which the electron motion in the transverse directions (y and 
z) is assumed to be completely frozen due to a strong confinement potential in both conduction and valence band(in 
other words, we consider only the lowest conduction and highest valence subband). Our results are likely to apply 
also to the case where the resonance occurs via some other intermediate state, but we have not considered this case 
explicitly. 

The longitudinal (x direction) motion is assumed to be free without defects or disorder. An important feature of 
quantum wires is the Coulomb interaction, which is typically unscreened and leads to a long ranged interaction with 
a well-known characteristic form 2 — involving the transverse wavefunction^S. The scale dependence of the unscreened 
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Coulomb interaction complicates considerably the analysis of Luttinger liquid formulae, but as will be seen also leads 
to the appearance of additional structures in the predicted RRS spectra. In this section we derive an expression for 
the RRS cross section of a general model with arbitrary electron interaction. Subsequent sections present results for 
short ranged interactions (where the analysis can be carried through in considerable analytical detail) and the more 
physically relevant long ranged interaction, which requires additional approximations. 

The appropriate theoretical starting point for analysis of Raman scattering is the following general Hamiltonian 

-fftot = Hy + Hll + H cv + Hjjf, (1) 

where Hy and Hll are the Hamiltonians for the valence and conduction bands respectively, and H, gives the 
coupling of these carriers to externally applied radiation. H cv describes the excitonic interaction between conduction 
band electrons and valence band holes, and may be important in certain conditions^!. We will assume excitonic effects 
are irrelevant in the energy regime of interest and will therefore set H cv = throughout this paper. In this case, 
the RRS process involves only one electron excited out of the valence band, so that interactions in the valence band 
can be also neglected. We consider a single one-dimensional conduction sub-band and a valence subband, linearize 
the dispersion about the Fermi level and include an (at this stage arbitrary) interaction between conduction band 
electrons. These considerations imply 

Hv — 2_j ^,p V r,p,s v r,p,s (2) 
r,p : s 

H LL =^V F (r P -k F )4 >p>s C r!p>s + — J2 y (9) C r 1 , Pl -g/2, Sl C L P2 +g/2, S2 C r 2 ,p 2 - 9 /2, S2 C r. 1 , P1+9 /2, S1 , (3) 

r,p,s T\ ,r 2 ,si ,s 2 q,pi ,f> 2 

where the conduction band Hamiltonian can be bosonized as in the standard Luttinger liquid theory^, yielding 

Hll = ( w p b l b p + w p°p°p) ■ 
p 

We approximate the valence band energy, E^ pl by a linear dispersion about the Fermi wavevector of electrons in the 
conduction band: 

E^ p w -O rrs - Vp(rp - k F ), 

where f2 rrs = Ep + Ep + E g is the RRS resonance energy (see Fig. [Ha)). Note that in this expression it is assumed 
that the valence band is also one dimensional. If transverse motion in the valence band is important, these degrees of 
freedom should be integrated out, which will broaden the valence band propagator. (Note that the structure of the 
quantum wire system means that momentum transverse to the wire need not be conserved in an optical transition). 
cl p s and v\ s are creation operators of electrons of chirality r = ±1 (left/right moving), wavevector p and spin s, 
in conduction band and valence band re spectively; anc j <jt are charge boson and spin boson creation operators in 
Luttinger liquid theory (see Refs. |l(Tl32| for a general review). vp(kp) is the Fermi velocity (wavevector) of conduction 
band electrons. V(q) is the effective ID electron-electron interaction within the conduction band. In Eq. |JSJ, the 
charge (spin) boson energy, u>P(u) p ), is related to the interaction V(q) via™ 



<Hp|Wi + — , (4) 

lj% = \p\v F , (5) 

where we have assumed that electron-electron interaction is spin-independent so that the spin boson velocity is the 
same as noninteracting electron Fermi velocity. 

Finally, for the electron-photon interaction Hamiltonia n, -Hmt , we consider only the three-leg vertex for which 
photon number is not conserved (see Fig. QJb) and Refs. [3fjl33ll34l ] ) . and neglect the four- leg vertex where photon 
number is conserved (see Fig. [He)), because the contribution of the latter does not give rise to resonance effects and 
is thus^S much smaller than the contribution of the three-leg vertex in near resonance conditions. We also represent 
the radiation by a classical field, so that 
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H int = (9lS,,,> +928 s -s>) {[4,p, S (*H,p-fc, S '(*) + V i,p,s>( t ) C r,p-k, S { t )V 

+ [4, P ,*(*H, P+M '(t) + < PjS '(*)cr, P +fc >s (t)]e iwt } , (6) 
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where u> and k are energy and wavevector of the photon interacting with electrons, gi and gi are coupling constants 
for non-spin-flip and spin-flip scattering respectively, and their actual values are not important in our study. In the 
remainder of this section, we present an explicit derivation of expressions for the non-spin-flip (polarized spectrum) 
RRS cross section following from Eqs. J2J, Q and 

Following our earlier work (Ref. 25.30|), we use second order time-dependent perturbation theory to calculate the 
rate at which H-^ causes transitions from the ground state |0) to some state |n) in which the valence band is filled 



int 

and the conduction band state has changed. One finds 



W{q,v,n)= Mm i£ [ T/2 d tl [ U to{ n \H^ v '\t 1 )H£ /2 ^' ,/a {t*)\0) 

n J-T/2 J-T/2 



(7) 



where q and v are photon wavevector and frequency shifts after Raman scattering and f2 is the mean frequency of 
incident and scattered photons during the process. We choose the backward scattering channel so that all wavevectors 
involved are along the wire direction for simplicity. Because the valence band is filled in the ground state, the part of 
the correlator involving valence electrons is simple: 

( W r 1 ,p 1 -g/2, Sl (*2) u r 2 , P2 , S2 (ti)4 3)J >3, S3 (*l)' ; r4,P4-g/2, S 4(*2))o 

= ^ 1 ,r 2 ^, r4 ^ 1 , P2+ ,/2^, W - 9 /2^ l!S2 ^, S4 e<« (ti " t;) e<-3( t 1 -^. (8) 

Using the space-time translational symmetry, Eq. (JJJ can be further simplified by representing fermion operators 
in coordinate space: 

pL />oo 

W(q, v; ft) = lim dR dTe l(lyT - <?i?) (O t (i?,T)6(0,0)) (9) 



(10) 



L — >oo 

Here (• • • )o are the expectation values on the ground state wavefunction and 



0{R,T)=y\ dx dt(f>{x,t)il; r . s (R + X /2,T + t/2)i>i s (R-x/2,T-t/2), 
„ . Jo Jo 



and 



0( X;t ) =i_J2e i ( E * t - px) = e^-^'+^^Six + rv^t). (11) 
p 

Eqs (f?Tf>— (| IL IL|) are our fundamental results. They show that the RRS process creates an electron- hole pair separated in 
space by x and in time by t, with the amplitude for a given space-time separation controlled by the function <j>(x,t), 
which is essentially the propagator for the valence-ban hole. Far from resonance (|fi — f2 rrs | ^> t/pfep), <j>(x,t) is 
short ranged in both x and t, so that O becomes similar to the ordinary density operator and W becomes a density- 
density correlation function 29-30 . As the mean photon energy is tuned closer to the resonance condition, (j> becomes 
longer ranged, so that O becomes nonlocal in both space and time. This non-locality will be seen to give rise to the 
interesting resonance effects, by allowing the light to couple to something other than the dynamical structure factor. 
We also observe that up to this point we have not made use of the one dimensionality in any important way: the 
equations may easily be generalized to two and three dimensions. Previous theoriesiS'iS'i^^^iSiiSSiSfiiSSiSSiS 3 . of the 
RRS process have been based on similar expressions but with a function cj> which essentially forces t = (i.e. no 
retardation effects). We shall see below that the time dependence is very important. 

We now incorporate the special features of one-dimensional physics by using the standard 3 ^ bosonization represen- 
tation of electron operators in Eq. (|10f> 

irk F x I rzr ( 

lMx,f) = ^=exp i V J— \e- pa / 2+trpx [cosh6 p b rp {t) - sinh9 p bl rp (t) + t 
\I2-Ka \ ~~i V P L L 




where exp(— 26 p ) — ujP/pvp = \Jl + V(p)/ttvf is the momentum dependent LL exponent, and a — > + is a conver- 
gence factor. Substituting Eq. (JllJ into Eq. I|1(J|) and using the following identity for linear boson operators, A and B 
(valid when [A, B] commutes with both A and B, the colons denote normal ordering and the <>o denotes expectation 
value with respect to the ground state of H LL ): 



e A e B = (e A e B ) : e A+B :, (12) 
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we can rewrite O as 



poo 

0(R,T) =LJ2 dte* n - nr "*Gr,.(-rv%t,t) : e*"" 



, p (R,-rv%t;T,t) .. e ® ra ,,,(R,-rv%t;T,t) 



where we have separately normal ordered the charge and spin boson operators and G rcr (x,t) = 
lim a ^o(27TQ!) _1 (ip r(T (x, i)V4 CT (0, 0)) is the electron Green's function in the conduction band. The phase operators, 
3>p and are respectively 

$,. s , p (i?, x- T, t) = 2 e ~° p/2 \Ft { - sinh 6 v sin\p(rx + v§t)/2] [bl rp e l ^ rR+ < T ^ + h.c] 



p>0 



+ cosh 9 p sm[p{rx - v p p t)/2] [bl p e-' lp{rR ~ v ^ + h.c] } , (13) 
$ rs , CT (R, x- T, t) = 2sJ2 e~ ap/2 ,£{sin[p( ra - v v F t)/2] [al p e'^ rR ~ vV ^ + h.c] } . (14) 

To calculate the expectation value of a product of normal ordcrings, we use two additional identities for the linear 
bosonic operators (valid under the same condition as Eq. (|12ll ) 

(eV) = (ene [AB1/2 (15) 

and 

(e A )o = e <A= )o/2) (16) 

so that 

(: e A :: e B :) = er^^e-^ 2 (e A e B ) 

= e -(A 2 ) /2 e -{B 2 ) /2^A+B^ oe [A,B]/2 

= e< AS >°. (17) 

Combining all the results from Eq. @ to Eq. (|17fl and defining the average photon frequency relative to the resonance 
frequency £l rrs as f2 = Q — £l rrs , we obtain 

ri.r-2 si,s 2 

/OO p pOC rOO 

dT dRe wT -^ R dhf^ih-n) dt 2 f r2S2 (t 2 ;n) 
-oo J JO JO 

x {exp [{$na uP (R, -riVph;T, ii)$ r2S2iP (0, -r 2 Vpt 2 ; 0, £ 2 )}o] 



where 



x exp[(* ri . ljCr (iJ,-nt;^ti;T,ti)# ra ^ i<r (0,-r2^t2;0 > *2))o] - 1} , (18) 



/„(t;n) = e l0t G rs (-r^i,t). (19) 



The — 1 in the last line of Eq. I|18|) arises from normal ordering^. 

The charge and spin part of the expectation values in Eq. 118(1 can be calculated respectively to be (at zero 
temperature) : 

(® risup (R, xi,T, h)$ r2S2iP (0, x 2 , 0, t 2 ))o 

= 2nr 2 ^ ^. e -^/2 e -i P v P T { sin[p(xi _ n vPtx)/2] sm[p{x 2 - r 2 v p t 2 ) /2]C rup C r2<p ^ R 
Jo P 

+ BintpCn + rii^*i)/2] sin[p(:E 2 + r 2 v p t 2 )/2]C^ ruP C. r ^ p e' lpR } , (20) 
where C r . p = cosh 9 p (sinh p ) for r = +1(— 1), and 

($ ri ai,ff »l! r ) ^l)*r 2S2 , o -(0, X 2 , 0, t 2 ))o 

= 2sis 2 5 ri , r2 /°° ^ e - Qp/2 sin[p(xi - w F ti)/2] sin[p(x 2 - v F t 2 ) /2]e ir '^ T (21) 

.A) p 

Eqs. (|18|I - H21|I represent a complete solution to the RRS cross section of a Luttinger liquid in the absence of excitonic 
effects. We will evaluate them analytically and numerically in the following sections. 
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III. LEADING ORDER EFFECTS IN BOSONIC EXPANSION 



A. Overview 



Although Eq. I|18(l is a complete closed-form expression for the RRS cross section, it requires a numerical evaluation 
which is not simple because final results are obtained by cancellation of rapidly oscillating terms. It is instructive to 
consider the analytical approximation obtained by expanding the exponential factors in Eq. (|18|l in a Taylor series in 
($$)" (n is an integer) and then evaluating the series term-by-term. Note that in the present problem (unlike in the 
evaluation of the electron Green function) each term in the expansion is finite because the sine function in Eqs. I|2()() 
and (|21[l removes the 1/p divergence in small p region and, by oscillating, ensures the convergence at large p. The 
expansion approach has a simple physical interpretation: the term of order n in the expansion corresponds to a final 
state with n excited bosons. 

Subsection B shows results obtained on the assumption of short-ranged interaction between electrons in conduction 
band, i.e. V(p) — Vq with an appropriate momentum cut-off. The well-known analytical methods of LL theory can 
then be applied, yielding a physical understanding of the features of RRS spectrum. Subsection C presents results 
obtained for the unscreened Coulomb interaction. 



B. Results for short-ranged interaction 



For simplicity here we assume that the valence band electron velocity at the Fermi momentum of the conduction 



band, 



F ' 



is much less than the conduction band velocity, vf, and can be neglected. In this approximation (which 



does not affect any essential results), 



„iQ.t 



f(t,n)=e lUt G rs (x = 0,t) 



where 



1(^ + ^.-2 

4 V vf v p 



sinh 8„ 



(22) 



(23) 



is the Luttinger exponent and Eq is an energy cut-off above which the Luttinger liquid model ceases to describe the 
physics. (Note that here we follow the standard convention in the literature^ and this a is nothing to do with the 
convergent factor used in Eq. Eq is expected to be roughly of the order of the conduction band Fermi energy, 

Ep^Zk. The time integrals can now be reduced to gamma functions, and the expansion evaluated order by order. 

We first show the results for the one-boson contribution in the non-spin-flip process. Here the spin boson term is 
canceled by the spin sum as required for the conservation of angular momentum. We have 



where the function A\ is 



W 1 (q,v;Sl) 



Ai(Ct,p, a) 



r 2 (- a ) 

2qvFV p 



Ai(n,g,a) e 26 »5{v - qv p ), 



' n-qv p /2 \ a _( + gv p /2 \ ' 
- { E 



En 



(24) 



(25) 



Thus the one boson term is a delta function at the plasmon energy, i.e. it gives the conventional CDE contribution to 
the scattering cross section. The amplitude for excitation depends on the difference of the mean photon frequency from 
the resonance condition and on the transferred momentum q. Far from resonance, cu q = qv p <C and W\ cx |f2| 2a ~ 2 , 
while at resonance Vt = and W\ cx sin 2 (7ra/2). Thus LL effects enter the CDE part iv ~ qv p ) of the spectrum in 
two ways (for short-ranged interaction): first, far from resonance, they changes the frequency dependence of spectral 
weights from (i~ 2 , the noninteracting result, to VL~ 2+2a . Second, near resonance (Cl < v p q/2) Luttinger Liquid effects 
resolve the singularities which are found in the standard Fermi-liquid expressions for the matrix elements: see Sec. [V] 
for a more detailed discussion. 
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For the second order (i.e. two-boson) contribution, we obtain 

T 2 (-a) f f dp 



W 2 (q,v,n) 



VF 



o p(q-p) 



A%(D,,p,q-p,a) cosh (20 p )8{v - qv p ) 



dp 



q p(p - q) 
q dp 



A%(£l,p,p- q,a) sinh 2 (29 p )5(i> - pv p - (p - q)v p ) 



o p{q-p) 



AZ(Q,,p,q-p,a) 8{v - qv F ) 



(26) 



with {a=p, a) 



^ - PlVg/2- P2Vg/2 
^ - PlVa/2 +P2Vg/2 




(27) 



Note that for large (l (far from resonance) W% ~ VL~ A+2a so that in this limit the two boson term is small compared 
to the one-boson term, confirming the validity of the expansion in the far-from resonance region. The first term in 
Eq. H26fl gives a renormalization of the strength of the 'CDE' pole. The second and third terms produce new affects 
appearing at second order: first near the CDE peak, branch mixing (second line of Eq. J2HJ0 processes appear, leading 
to a continuum absorption beginning at the CDE threshold, v = qv p . Secondly the spin singlet combination of spin 
bosons (third line) can be excited via the two spinon, (aa^aa) process, giving rise to the so-called "SPE" mode at 
v = qvF- 

In Fig. we show the perturbatively calculated LL polarized RRS spectra including one and two boson contri- 
butions to the Raman intensity, for a particular choice of Luttineger exponent The data are plotted as a function of 
energy transferred at fixed momentum transfer q; the different graphs show results for different values of the laser 
frequency relative to the resonance condition. One observes that (i) the overall spectral weight decays very fast as 
51 is moved off resonance, and (ii) the "SPE" peak is only noticeable near resonance, (iii) The continuum structure 
arising from the second term in Eq. (|26[) is visible but not very distinct. We also remark that the three boson term 
in the expansion (not calculated in paper) will produce an additional continuum structure lying between the "SPE" 
(lu = qvp) and the plasmon (u) — qv p ) energies, due to the coupling between charge and spin bosons via the RRS 
process. 



C. Long-ranged Coulomb interaction 

We now consider the modifications arising in the physically relevant case of a long range, unscreened Coulomb 
interaction. As discussed elsewhere^ this corresponds to a scale dependent Luttinger interaction parameter, a q = 
ao In 1 / 2 (qo/q) (for q -€1 qo) and the charge mode plasmon energy^, 

LU p q ~4a VFqln 1/2 (q /q) (28) 

Here qo ~ 2.5 /a (a is the typical width of QWR) is the momentum scale below which Coulomb effects become 
important and ao = y/ e 2 /2eow_F7r^i for typical QWR structure. In the GaAs QWR material which was recently 
studied^, typical LL parameters are ao ~ 0.389 and qo ~ 1.4 x 10 6 cm -1 . 

The scale dependence of the interaction prevents analytic evaluation of the Luttinger liquid formulae, but the 
bosonic expansion can still be carried out, although both the integrals defining f(t, Cl) and those defining the <& must 
be evaluated numerically. 

In Fig. |3we show the one and two boson contributions to the polarized RRS intensity of ID QWR using Coulomb 
interaction parameters corresponding to the system studied in Ref^. The basic features of the results are similar 
to the short-ranged interaction, having one weak singlet spinon peak at the "SPE" spinon energy u> = qv a — qvF 
arising from the two coupled spin bosons, and one strong charge boson peak at the plasmon energy cu = qv q . Their 
relative strength varies according to the resonance condition. Off resonance, the two-boson-singlet-spinon peak is 
much weaker than the charge boson plasmon peak, while their total spectral weight is also very small compared with 
the near resonance result. 
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Near resonance, an additional feature appears: the two-charge-boson contribution gives another resonance peak 
at uj = %Uq/ 2 > w Ji arising from the nonlinear q dependence of io p q due to the long-ranged Coulomb interaction. 
The resulting curvature acts to broaden the spectral weight on the higher energy side. Next order bosonic terms 
will contribute the spin-charge mixing weights from {oop, poo) type higher order correlations. This will introduce 
a continuum of relatively weak spectral weight which is located in the energy region between the spinon and the 
plasmon energies, as noted previously. 



IV. BEYOND THE BOSONIC EXPANSION 



We now consider effects beyon the bosonic expansion in the analytically tractable finite range interaction model. 
We begin with the exact expression, Eqs. i|20|l - i|21|) . The values of mean position R and time T are expected to be 
roughly inversely proportional to the momentum q and frequency v transferred to the system. If v and q are small 
compared to the basic scales {Ef and Uf) we may use a long wavelength approximation to evaluate the phase factors 
in the standard way, obtaining 

/>oo />oo />oo />oo 

W^g( q ,v,n)= dhf; isl (h;Q) dt 2 f r2S2 (t 2 ;Cl) dTe™ T dRe~^ R 

JO JO J-oo J-oo 

x { [F^ ri {r x R - v a T- t u t 2 T° [F r p ^ r2 (R - v p T: h,t 2 )) a ^ 

x [F? u r 2 {-{R + v p T)-MM)f? -l}, (29) 
where the exponents are, a a = ^5 ri} r 2 sis 2 , or p = ^C ri C r2 , and offi = iC_ ri C_ r2 . The F function is: 

F^XMM) - {v a itl _ v a 2t2) 2_ [X + lar (30) 

where = (v a +rvp)/2 for a — o, p (v a — vp)- 

We have not attempted a direct evaluation of this expression. Instead, we note that the spectrum is expected 
on general grounds to consist of two delta functions, at the CDE and SPE energies, and two continua. The delta 
functions are the features most easily observable experimentally, and it is possible to obtain convenient expressions 
for their weights. We note that terms giving rise to delta functions at the CDE energy or SPE energy must be 
functions only of R — v p T and R — v a T respectively. To isolate these terms we apply the identity ABC — 1 = 
(A - l)(B - 1)(C - 1) + (A - 1)(B - 1) + (A - \){C — 1) + (B — l)(C - 1) + (A - 1) + (B - 1) + (C - I) to Eq. 
(|29J) . Only terms involving just one of the three F- functions can contribute to the delta function, so we find that for 
v. q > the intensities may be written 



Icde{&) =W p (n)5(u-v p q) + ... (31) 

(32) 



I S p E (n)=l J2 s lS2 W^(Cl) j 5{u - v a q) + ... 

\S1S2 = ±1 / 

where the ellipses denote the continuum terms and the delta function coefficients are, respectively: 

/>oo />oo />oo 

W p (q,Q) = 8n dhf^h) dt 2 ft 2S2 (t 2 ) dR'e^ qR ' {[F^R 1 - t u i 2 )p - 1} (33) 

Jo Jo J — oo 



W^-(q,n) = 2tt / dhf^h) / dt 2 ft 2S2 (t 2 ) / dR'e~^ Rl {[F^iR-MM)]^ - l} • (34) 

JO Jo J-oo L ' 

These formulae have an involved analytic structure, which includes a discontinuity when either the incoming or 
the outgoing laser frequency is precisely in resonance with the valence band to conduction band Fermi level energy 
difference. For large frequencies, the correct asymptotic behavior arises from cancellation among oscillating terms. To 
make the structure manifest and obtain forms which are convenient for numerical evaluation we employ the further 
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series of transformations given in the Appendix, leading for > v p ^ a qj2 to 

2l ' " ' 



W p (q, n > v p q/2) = 
(( 



Qsin (£(i+£))r(l-2«) r /4^ r /a^ cos 2g (cotg) ^ 



ft 1 " 2 " 



v~ 2 " / 



COS( 



v P q I 1+cos20cos26> 
20 V 2 



VI + cos 20 cos 29 
1 



cos < 



v r hoq / 1+cos 20 cos 20 

"iov 2 



(35) 



(( 



22+«r(l - 2a) 
fii~ 2 « 

1 



tt/4 



ir/2 



/0 

l-2a 



COS — 



v&q 1 1+cos 20 cos 20 



2<> 



d<9- 

V 



cos26> 



y/1 + COS 20 COS 26* 

1 



,- 2a 



COS ( 



v^q 1 1+cos 20 cos 29 
2STV 2 



For 



< Vp^q/2 we obtain 







2 f ~ Q / 1 + n . r /4 1± r' 2 ,„ cos 2 0(cot0) a 

< v p q/2) = sin ( n—^— ) T(l - 2a) / <20 / dO 



( \l-2a 

{v P q) 



o VI ~~ cos 20 cos 20 



1 



l-cos20cos 28 | 2S1 , 
2 + u p <j 1 



1— 2a 



l-cos20cos26> _2f^ , 
2 -Up<j 1 



l-2a 



and 



7r/2 



cos26> 



Ii; r (<y. 5! < ^g/2) = 2 " ° r(l -2a) \ ' d 

(Ufft?) Jo Jo V 1 - cos 2 cos 29 

( 



l-cos20cos 26 , 2Q. 

2 + u„g 



sine 



l-2a 



l-cos20cos 2fl _ 2Q 
2 u CT q 



S11K 



l-2a 



/ 



The change in analytic structure occurring at 



(36) 



(37) 



(38) 



tor of v 7 ! i cos 20 cos 20 in the denominator. In the 

, 1-2q 

l-cos20cos 2fl _j_ 2Q 



Vp^q/2 is manifest in these expressions by the fac- 
Q 



2 'Up^g 

singularity at = = 0, which is absent for 



< Vp.aq/2 regime ^/l — cos 20 cos 29 and the 
have singularities at the same point, leading to a strong (but still integrable) 



n 



> v PtCr q/2. To understand the singularity one may expand for small 



, 9, finding the leading behavior (see the Appendix) 

II", 



1 2Q 

P ,a^-J a ( ) ( <v p , a q/2) 

a v p , a q 



(39) 



with J a a function with an additional integrable divergence at fl = Vp. a q/2. The extra contribution arising from 

Eq. (|39|l leads to steps in the CDE and SPE intensities as the frequency is moved across resonance, with magnitude 
diverging in the non-interacting limit a — > 0. 

Fig. 0] shows, on a logarithmic scale the CDE and SPE results calculated from Eqs. (|35[l - i|38|) . For clarity of 
presentation we have removed the prefactors ((r(l — 2a)/ (v Pta q) a and have expressed the dependence on laser 
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frequency via the ratio of the laser frequency measured from resonance Q to the mode frequency v (note that at fixed 
momentum transfer q the SPE and CDE mode frequencies will differ). Far from resonance, the CDE intensity is seen 
to be much larger than the SPE intensity, but close to resonance, the two are comparable. The main panel of Fig. 
shows the ratio of the SPE to the CDE absorption. We see that the ratio is only appreciable if the incident photon 
energy is essentially 'on resonance'. The ratio exhibits a derivative discontinuity (seen more clearly in the inset) at 
|0| = v /2 (here v — Vp^q represents the resonance energy for charge and spin modes respectively); the on-resonance 
ratio is about 0.5 independent of the Luttinger exponent. (Of course the factors of velocity which have been removed 
from the results will lead to an additional dependence on a). 



V. RPA CALCULATION 



In order to compare our calculation to previous literature we present the analogue, for the model considered here, of 
the " Fermi Liquid' calculation discussed extensively in previous workiSjiSiiLiSiSiiSi. The previous calculations neglect 
backscattering entirely, and treat the forward scattering part of the interaction in the RPA approximation; we make 
the same approximations here. We obtain analytical expressions apparently not given in previous literature. We 
consider here a model with short ranged interaction, parametrized by a constant amplitude V. The CDE feature is 
then a zero-sound collective mode with a velocity shifted from the Fermi velocity by the interaction strength. 

We evaluate the diagrams shown in Figl^l The analytic expression corresponding to the diagram labelled 1 is 

D1 = -TJ2 [ Y G ^ m + iv / 2 + iLJn ' p + ll 2 ) G ^P + iu » + iv ) G ^ *fn) (40) 

Setting the valence band velocity to zero, analytically continuing on the laser frequency £1 and measuring it from 
resonance and performing the frequency sum yields 

( _ 

fp 1 fp+q 1 



D\ 



dp 1 

2ir ii/ + e p - e p+q 



+ if/2 + e p ) (Jl - iu/2 + s p+q ^ 



(41) 



where f p is Fermi distribution function. We next linearize the conduction band dispersion, perform the p integral 
(bearing in mind that we must take the principal value), drop terms of order vq/Ep and analytically continue on 
v. Unlike the Luttinger liquid expressions obtained in previous sections, this expression as it stands is infinite in 
the range ft < vq/2. Previous work2&Si resolved this divergence by introducing a phenomenological broadening 
parametrized by a quantity A. We follow this procedure here, but emphasize that because the behavior in the region 
51 < v /2 is determined entirely by this phenomenological parameter, our results in this region have, strictly speaking, 
no meaning. We obtain 

'taiW^Utan- 
Dl = n>, q) ( ^ ^— -L I , 42 ) 




Here 



n(^,g) = ^— f - s + 7^]=— — 5 1 o s (43) 

Zirvp yv — VFq — io v + VFq — ioJ irvp v*- — Vpq z 



D2 = (D22) 2 — - — (44) 



with 



The analytical expression corresponding to the diagram labelled 2 is 

V 

i + xv 

D22 = -T^2 I Y G d( in + il/ / 2 + iuJ n,P + q/2)G c {p + q,iuj n + iv)G c (p,iuj n ) (45) 
-^G c {p + q, iuj n + iu)G c (p, iu> n ) 

(46) 



{vFqY 



irvp v 2 — vpq 2 
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The arguments leading to Eq. I|42|l may be repeated here. The divergence is weaker; indeed the obtained expressions 



are finite everywhere except at 
as above yields 



v/2, where there is a logarithmic divergence. Introducing the same broadening 



D22 = 



dp 



1 



fp+q — 1 



2tt iv + e„ — £ 



(fl-ii// 



2 + e 



p+q , 



= n(g, I /)-ln 



The total unpolarized RRS cross section, in this approximation, is thus 



A 2 



SI -v/2) +A 2 



I 



RRS 



U(q, v) 



tan 



-i ( n+v/2 
A 



tan 



1 / Q-v/2 



ttvf v 2 — v'pq 2 



tan 



A 



i / n+i//2 




lln 2 



(n+y/2) 2 +A 2 
(n-y/2) ; "+A 2 



— tan 



i f n-u/2 
A 



A 




1 + 



tt | In 2 



/ 

(0+i//2) 2 +A 2 
(f2-i//2) 2 +A 2 



IXVp 



v 2 — v 2 q 2 



(47) 



(48) 



with Up = Vp y 1 + Observe that if the phenomenological broadening parameter A is set to the RPA approxi- 
mation to the RRS intensity diverges when O = ±u/2, in other words when the incoming or the outgoing photon is 
in resonance with a transition of the system. 

From these formulae the SPE and CDE intensities may easily be obtained. We have 



Icde 



lspe 



1 



hi" 



(fl + v/2^ 2 + A 2 
fl - v/2\ 2 + A 2 



5(v-qv p ) 



2v F 



tan 



-1 ^l+v/2 \ 



tan 



_X ( n-v/2 
A 



A 




+ +A 

- \2 
n-v/2) +A 



b(y — qVF) 



(49) 



(50) 



The CDE/SPE ratio is plotted in the inset to Fig. [5]for A = O.Iia We see from this inset that the ratio is generically 
very small (smaller than the corresponding ratio in a Luttinger liquid), but that for all Q < u p . CT g/2, the divergence 
of the SPE term (cut-off here by a small value of the phenomenological parameter A), yields a large value of the ratio. 



VI. DISCUSSION 
A. Overview of calculation 



In this paper we have presented a complete analytical theory of Resonant Raman Scattering in Luttinger liquids, 
and for comparison we have also given the analogous results for what is referred to in the literature as the 'Fermi 
liquid' approximation. In this section we discuss the implications of our results. In this first subsection we outline the 
essential assumptions underlying the calculation. In the second subsection we discuss the relevant parameter regimes 
and the different behavior expected in each regime, and in the third subsection we discuss the extent to which Raman 
scattering has observed, or can be expected to observe, evidence for Luttinger liquid behavior in quasi one dimensional 
systems. 

We note that both Luttinger liquid and 'Fermi Liquid' calculations require three assumptions: 
(i) that the energy v and momentum q transferred to the one dimensional electron system are small in comparison 
the the Fermi energy and Fermi momentum of the one dimensional electron system 
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(ii) that the Raman process involves excitation from a one dimensional valence band to an empty state in a one 
dimensional conduction band, followed by decay of an electron from a filled state in the conduction band into the hole 
state in the valence band 

(Hi) that 'excitonic' correlations between conduction band states and the intermediate state may be neglected. 

Assumption (i) is essentially a condition that the experimental resolution is sufficient to reveal the low energy 
physics of the system of interest. From the theoretical point of view it could easily be relaxed at the expense of 
introducing a more complicated description of the conduction band. 

Assumption (ii) is clearly applicable to strictly one dimensional systems such as carbon nanotubes (where the valence 
band is clearly one dimensional) but may not be applicable to quantum wire structures created by lithographic or 
MBE techniques on a three dimensional substrate. In this situation, valence band carriers may be able to move 
transverse to the wire. Because momentum transverse to the wire is not conserved in the optical absorbtion process, 
in this case the function <f)(x, t) in Eq. Ijllfl must be integrated over a range of transverse momenta, leading to a 
broadening of the resonance. However, assumption (ii) may also be relaxed. The key issue is that the intermediate 
state (of whatever origin) disperses only along the wire. 

The crucial assumption is (Hi), neglect of excitonic correlations. These are likely to play a crucial role in the 
QWR structures, keeping the valence hole near the wire and therefore allowing a sharp resonance to occur. Neglect 
of excitonic correlations entered the theory at several points, and the theory cannot be trivially modified to include 
them. Construction of a theory including excitonic effects is an important open issue in the field. 

By combining the three assumptions given above with standard second order perturbation methods, we obtained 
results, presented in Sec. IIVI and Eqs. (|35|) - (|38[) . for the resonant Raman scattering intensities of a one dimensional 
electronic system described by the Luttinger liquid model. For comparison, we also presented analytical expressions for 
the Raman intensities of the previously discussed 'Fermi liquid model" (in which the only electron-electron interaction 
considered is forward scattering in the density channel, and this is treated in the random phase approximation (RPA)). 
In the far-from resonance regime, we showed that a perturbative treatment could be applied, and this was used to 
provide a detailed characterization of the different features of the spectrum, including 'CDE' and 'SPE' peaks and a 
continuum absorption and the dependence of the intensity of these features on the difference of the laser frequency 
from resonance. We also obtained simple, easily evaluated expressions for the integrated intensities in the SPE and 
CDE peaks for all values of the laser frequency, and from these we obtained results for the ratio of SPE to CDE 
intensities. The expressions have an interesting analytical structure in the near-resonance region, which is discussed 
below. 



B. Physics of RRS spectrum— difference between Raman spectra of Luttinger and Fermi Liquids 

The Raman process creates a particle and a hole in the conduction band, at a space-time separation controlled by 
the difference of the laser frequency from a resonance condition. The particle-hole pair is in general not an eigenstate of 
the conduction band Hamiltonian, but decays into eigenstates. The important question, therefore, is "what properties 
of the eigenstates are reveals by Raman scattering" . In the one dimensional context one may sharpen this question to 
" what aspects of the Raman spectrum reveal characteristic features of the Luttinger liquid physics expected to occur 
in one dimensional systems?" . In order to answer this question precisely, it is necessary to discuss the eigenstates of 
one dimensional electronic systems. 

We are concerned here with polarized channel Raman scattering in systems with negligible spin-orbit coupling 
and therefore must consider excitations which do not change the total spin of the one dimensional system. In a 
Fermi liquid, the eigenstates are particle-hole pairs, cither organized into collective models or existing as incoherent 
" particle- hole continuum' states. Further, in a Fermi liquid at low energies and long wavelengths one may restrict 
attention to states consisting of a single excited (quasi)-particle or (quasi)-hole: multipair states are important only as 
virtual processes renormalizing matrix elements and dispersions. In a Fermi liquid with repulsive interactions, there 
are two kinds of relevant states. One is the zero sound (or 'plasmon') mode, referred to in the Raman literature as the 
'CDE' mode. (In principal other collective modes exist; these are rarely important in practice). The CDE excitation 
typically has a well defined energy- momentum dispersion relation v — v p (q) \ q\ with velocity v p (q) greater than the 
Fermi velocity Vp. 

The other class of states, referred to in the Raman literature as 'SPE', are the particle-hole continuum states. In two 
and three dimensional systems, these excitations exist in the range < v < VFq', in d = 1 they exist only in the range 

2 2 

v f\i\~ — \ v \ — v f M + if^ - with p — d i e ^g^i — a measure of the curvature of the quasi-particle dispersion. 
In any dimension, the 'SPE' or continuum excitations make an important contribution to the specific heat and to most 
response functions. However, the contribution of the continuum excitations to the structure factor (density response 
function) is typically much smaller than that of the CDE (plasmon) mode. In d > 1 the CDE contribution to the 
structure factor is larger than the SPE contribution by a factor of the dimcnsionless long wavelength interaction 
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strength, which for a charged Fermi liquid involves factors of (<7tf/<?) ~ arising from the long ranged nature of 
the Coulomb interaction (here qtf is the Thomas- Fermi screening length). In d = 1 the factor arising from the 
Coulomb interaction is only logarithmic; however the special kinematics of one dimensional systems lead to additional 
constraints. Within the RPA approximation the SPE contribution to the structure factor is smaller than the CDE 
contribution by a factor of (q/po) 2 . In the linearized dispersion limit (po — > oo) the SPE contribution to the d = 1 
structure factor vanishes entirely within RPA approximation, and indeed this vanishing may be shown via a Ward 
Identity? 8 to occur to all orders in the interaction. We emphasize, however, that in the Fermi liquid approximation 
to the one dimensional electron gas, the SPE excitations exist, and may be revealed either by considering response 
functions other than the structure factor, or from the specific heat. A specific example of a response function other 
than the structure factor is provided by the Raman intensity near resonance. 

In one dimensional systems, the Fermi liquid approximation is believed to provide a very poor description of the low 
energy physics, which is believed to be more accurately represented by the Luttinger liquid model. Consider first a 
spinless electron gas. This could be realized in practice by applying a magnetic field large enough to fully spin-polarize 
the conduction band. It is believed that a correct treatment (for example, via the Luttinger model) would predict 
that the only elementary excitation is a boson, roughly equivalent to the zero sound or plasmon mode. Thus in the 
spinless Luttinger liquid case, in all response functions and also in the specific heat, only superpositions of plasmons 
would be observed, whereas the 'fermi liquid' approximation predicts one would observe features at both energies. In 
Raman language, in the Luttinger liquid case only the CDE mode would be visible, no matter how closely the system 
is tuned to resonance, whereas for a hypothetical ID Fermi liquid, excitation at the SPE energy would become visible 
as the laser frequency is tuned to resonance. 

Now consider the Luttinger liquid formed by the usual electron gas, in zero magnetic field. In this case, two classes 
of elementary excitation exist: charge and spin bosons. States of total spin zero may be constructed from states of 
two or more spin bosons, even if no charge bosons are present, and these may make important contributions to many 
response functions and to the specific heat. In the linearized dispersion limit, states involving only spin bosons do 
not appear in the structure function, which couples only to charge. However, curvature in the dispersion leads to 
charge-spin coupling and in particular to terms allowing a virtual charge boson to decay into two spin bosons. Thus, 
just as in the Fermi liquid case, if nonlincarity in the underlying dispersion is neglected, only the CDE mode is visible 
in the structure factor, but at order (q/po) 2 SPE contributions appear, and of course SPE contributions appear in 
response functions other than the structure factor. 

To summarize this subsection: as long as the electron number is not changed, the apparently profound differences 
between a Luttinger and a Fermi liquid produce only minor, quantitative differences in the excitation spectrum of 
an electron gas in zero magnetic field. In both cases, one has two classes of excitation, which may be labelled CDE 
and SPE. Far from resonance, the matrix element determining the Raman spectrum is the density operator p q v plus 
additional operators which vanish as the ratio vjQ. vanishes (i.e. as the laser frequency is moved arbitrarily far from 
the resonance condition). In the limit v/fl — > the Raman spectrum is the same as the structure factor. In both 
Luttinger and Fermi liquids the SPE mode makes an 'intrinsic' contribution to the structure factor which is smaller 
by a factor of (q/po) 2 independent of difference of laser frequency from resonance and which vanishes when curvature 
is neglected. However, for a spin polarized electron gas, there is an important difference between a Luttinger and a 
Fermi liquid: in the Fermi liquid approximation an SPE branch of excitations exists, whereas in a Luttinger liquid it 
does not. 



C. Approach to resonance 

The previous subsections have shown that in the polarized channel of RRS experiment (i.e. no spin-flipping of final 
electron configuration), the differences between the predicted Raman spectra of a Luttinger liquid and Fermi liquid 
approximations to the one dimensional electron gas are quantitative, not qualitative. Very far from resonance, the 
Raman spectrum is determined by the structure factor, so differences between the two cases involve differences in the 
structure factor. These are unimportant, arising from different coefficients of the (q/po) 2 terms in the structure factor. 
Because the coefficients involve non-universal factors, this difference does not provide a useful distinction between 
the two cases. A more useful distinction arises from consideration of the changes arising in the spectrum as the laser 
frequency is brought closer to a resonance condition. As this occurs, operators other than the structure factor begin 
to contribute to the Raman spectrum, and the nature and coefficients of these operators may be used to distinguish 
between the two cases. 

In discussing the near-resonance behavior it is helpful to refer to Fig. [f)| which shows the diagrams commonly 
considered in the "Fermi liquid" treatment of the RRS process. From these diagrams one sees that in the 'Fermi 
liquid' approximation, the states created by the Raman process have nonvanishing overlap with the two sorts of exact 
eigenstates (CDE and SPE) of the system. Diagram 1 gives the dominant contribution to the probability of creating 
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an SPE excitation. The direct overlap between the state created by the Raman process and the exact eigenstate of 
the system means that in the absence of 'excitonic' correlations between the particle-hole pair and the intermediate 
resonant state (sketched as dash-dot lines in diagram 1 of FigEJ, this diagram diverges strongly. (The diagram might 
also diverge strongly even in the presence of excitonic effects: this point has not been investigated.) Indeed one sees 

from Eq. I|42[l that very near resonance the term diverges as (^l ± v/2^j and (in the absence of excitonic correlations 

or of the phcnomcnological broadening A) is infinite for Q < v/2. On the other hand, one sees from diagram 2 that 

the coupling of the Raman process to the CDE excitation goes through a range of virtual states (represented by the 
triangles in diagram 2); this broadens the resonance effect so that the divergence as resonance is approached is only 



logarithmic and the result is finite for 



< v/2. However, the one dimensional kinematics lead for O > v/2 to a 



much greater CDE amplitude than an SPE amplitude, so that the SPE/CDE ratio only becomes of order unity very 

close to the resonance condition. Further, for VL < v/2 the Fermi liquid approximation to the SPE intensity is, strictly 

speaking, infinite, and must be regularized by consideration of processes so far omitted from the approximation. At 
present there is no theory of this regularization, which is parametrized by a phenomenological broadening A. The 
'Fermi liquid' results presented in this paper in the near-resonance regime are probably meaningless. We suspect that 
this issue also arises in the higher dimensional calculations performed so far. 

In the Luttinger liquid approximation, the particle-hole pair created in the Raman process has vanishing overlap 
with any exact eigenstate. This has three consequences: first, a continuum absorption arising from multiboson 
excitations is predicted to lie in between the SPE and CDE peaks (actually, this continuum would exist also in a 
Fermi liquid, but would be much weaker). Also, away from resonance, the SPE/CDE ratio is larger than in the 
Fermi liquid approximation, and it varies with a power related to the Luttinger liquid exponent a. Finally, in the 

£1 < v/2 both SPE and CDE contributions are predicted to be finite, although both exhibit a 
= v/2. Fermi liquid behavior is recovered because the magnitude of the discontinuity is a -1 



on-resonance regime 
discontinuity when 



n 



and diverges in the noninteracting (a — * 0) limit. 



We suggest that in the on-resonance regime 







v/2 neither the Fermi liquid nor the Luttinger liquid calculation 

are likely to be quantitatively reliable, because excitonic correlations (neglected here) are likely to be important: near 
resonance the intermediate state lasts so long that it must interact with the conduction band excitations, and these 



interactions are likely to have a nontrivial interplay with the non-analyticities arising when 



v/2. 



VII. SUMMARY 



In this paper we have used the full Luttinger liquid model to analytically and numerically calculate the RRS spec- 
trum for both polarized and depolarized spectroscopy, and have presented for comparison a 'Fermi liquid' calculation. 
We obtained results for both a short ranged (screened Coulomb) interactions and for the more physically relevant case 
of the long ranged Coulomb interaction. We clarified the difference between the 'Fermi-liquid' and Luttinger liquid 
results, and argued that proper treatment of excitonic interactions (neglected in all treatments so far) is essential 
for obtaining reliable results. Our results in Fig. and Fig. shows that the RRS spectra of both short-ranged 
and long-ranged interactions are qualitatively similar, except that in the unscreened Coulomb case the momentum 
dependent charge velocity makes the phase space of multiboson excitations to be highly restricted, which shift the 
multiboson excitations to higher energy than the plasmon energy uj — qv q , causing such unusual broadening. The 
SPE-CDE ratio is generically larger in a Luttinger liquid than in a Fermi liquid, except on-resonance, where excitonic 
effects probably invalidate either calculation. 
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IX. APPENDIX A: EVALUATION OF DELTA FUNCTION CONTRIBUTIONS 

A. Overview 

The Appendix presents the additional transformations needed to put Eqs. (|33|l - l|34(l into forms convenient for 
further discussion and numerical evaluation. We have to deal with an expression of the form (note that for simplicity, 
we define Q to be the light frequency with respect to the resonance energy throughout this Appendix.) 

f°° dt 1 e^ n+l ^ tl f°° dt 2 e^ n ~ l£ ^ , 

W= ^-r+s / -*-i+5 I(ti,h) (51) 

Jo t 1 Jo t 2 

with e — > + being an infinitesmall converging factor and 

/oo 
dR e~ zqR F{t 1 ,t 2 ,R) (52) 
-OO 

and F given by 

where for the charge case v — v p and f3 = (1 + a) /2, while in the spin case, v = v a = Vp and (3 = 1/2. As in Sec. 
IIHI we have assumed that the valence band velocity is much smaller than the Fermi velocity in the conduction band. 
The analytic structure of the F function means that we can write (note we have also rescaled R) 

I(h,t 2 ) = 2sin(vr/3) f dR sin(^)F 2 (B, A, R) (54) 
J A 2 



with 



/ B 2 — R 2 \ ^ 

mi^R) = [^ w ^) (55) 

and B = t\ + t 2 and A = \ti — t 2 \. 

By making the following changes of variables: R — > \/u] u — > v + A 2 ; v — ► \/ B 2 — A 2 sin# and restoring the explicit 
forms of B and A we obtain 

, r' 2 (cos 0) 1+2/3 (sin 6) 1-2/3 sin ^-Ji:{ + f2 - 2ti*a cos 20 
I(t 1 ,t 2 ) = 2sm(irf3)t 1 t 2 dO y - ^ ^ /, ; , 2 V ' =1 — (56) 



■s/t[+q - 2ht 2 cos 20 



B. Evaluation, large O 



If O > -# then it is convenient to rotate the t integrals via 



getting 



with 



h -> i Tl (57) 
t 2 -» -ir 2 (58) 



W/ = / 1+5—/ i+s-^(n,r 2 ) (59) 



I(n,T 2 ) = 2sin(7r/3)riT 2 / d0 2 (60) 

JO V T 1 + T 2 + 2r l T 2 COS 261 
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It is then useful to define 



T\ = T COS ( 

T2 = r sin (f 



so that 



tt/2 ,ir/2 yoo dr 



al . ,s cos 9 1+2/3 (sin 60 ^sinh^ VI + sin 20 cos 20 
v ; '- '- '- - 2a VI + sin 20 cos 20 



(61) 
(62) 



(63) 



The r integral may now be done. It is convenient to shift the origin of the integral by 7r/4 and symmetrize in 
and to combine the two terms arising in the spin term, getting 



n 



cos < 



1 

v p q 1 1+cos 20 cos 26 



n JO 
\-1a 



2i+ Q sinf^2^))r(l-2a) r /4 r /2 

) 

\ 
/ 



VI + cos 20 cos 20 



(64) 



2!> 



2 



t- 2Q 



COS( 



-» p g / 1+cos 20 cos 28 



2H 



VI + cos 20 cos 20 



(65) 



1 



VfQ I 1+cos 20 cos 29 



2Q 



1 



/ 1+cos 2j^ cos 29 \ 
2Q V 2 I J 



) 



2 ')) 

This formula has integrable singularities at = = 7r/2 and 6 = tt/2 and is convenient for numerical evaluation. 



C. Evaluation, small Q 



If O < then the rotation cannot be made and we should instead define 

t% = t cos 
t% = t sin 

Transformations similar to those leading to Ea 1631 then yield 



(66) 
(67) 



and 



W P (Q) 



2i- a 



1 + a 



sin 7T-— r(l - 2a) 



7r/4 i-tv/2 



Jo 



d6- 



cos 2 0(cot0) a 
VI — cos 20 cos 20 



l-cos20cos20 i 2i1 



v p q 



sm < 



l-2a 



l-cos20cos2g _ 20, 
2 v p q 



S11K 



l-2a 



2§ r 2 r( i-2a) r% r /2 - c ° s2g 

(^q) 1 - 2 " V Vo Jo VI -cos 20 cos 20 



l-cos20cos 26 , 2Q 
2 + B„g ' 



sm < 



l-2a 



l-cos20cos 26 _ 2Q 
2 v„q 



S11K 



l-2a 



(68) 



(69) 
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We observe that unlike Eqs. (|64|l and (|65|l these expressions have a double singularity at small (j>, 9 and indeed lead 
to RRS spectra of quite different magnitudes. To see this point more clearly, consider the integrand, of Eq. 168|) 
as </>, 6 — > 0. We obtain 



^V¥TW y^^, + ^y 

Changing variables to 9 = C, cos ip and = C sin ^ we have 



J d^6d8I ch 



CdCdip- 



Integration over £ gives an answer proportional to l/a. 



l-2a 



/ \ l-2a 



(70) 



(71) 
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FIG. 1: (a) Schematic representation of RRS process in the direct gap two band model of electron doped GaAs nanostructures. 
LOi and ujf are the initial and final frequencies of the external photons, (b) and (c) are the three-leg and four-leg scattering 
vertices for electron-photon interaction. 
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FIG. 2: Polarized RRS spectra calculated via the bosonic expansion method for various resonance conditions, Q = Q — £l rrs . 
One- and two-boson contributions have been plotted separately in order to show their relative contributions (see text). Finite 
broadening factor 7 is introduced to express the delta function. Note that the overall spectral weights decreases dramatically 
off-resonance, as indicated by the individual scale factors on the right hand side of each plot. 
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FIG. 3: Polarized RRS spectra calculated via bosonic expansion method for various resonance conditions in the LL model 
with long-ranged Coulomb interaction. One- and two-boson contributions have been plotted separately in order to show their 
relative contributions (see the text). Finite broadening factor 7 is introduced to plot the delta function contribution. 
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FIG. 4: Logarithm of polarized-channel Raman scattering cross section (normalized to appropriate mode velocity), Irrs, 
plotted against incident laser frequency (normalized to mode excitation frequency v — v PtC ,q and measured from the average 
of incident photon resonance energy and outgoing photon resonance energy with respect to the resonance energy, Q rrs ), for 
two different values of the Luttinger exponent a. The plotted intensities are coefficients of delta functions describing coherent 
charge (CDE) and spin (SPE) final states and are computed via Luttinger liquid methods as described in the text 
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FIG. 5: Comparison of ratio of delta coefficients for CDE and SPE absorption for Luttinger (solid and long-dahsed lines) and 
RPA-Fermi liquid (dash-dot and short-dashed lines) model, normalized to appropriate mode frequencies and plotted against 
incident laser frequency (scale as in Fig. Inset: same ratio, displayed in linear scale. 
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FIG. 6: Feynman diagrams considered in evaluation of 'RPA' approximation to unpolarized RRS cross section. Diagram 1 gives 
the dominant contribution to the SPE absorbtion. Excitonic interactions (such as those shown as dash-dot lines in the lower left 
corner of this diagram) are not considered in this paper. Diagram 2 gives the dominant contribution to the CDE absorbtion. 



